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Abstract 

We have searched for neutrinoless r lepton decays into three charged leptons using 
an 87.1 fb^^ data sample collected with the Belle detector at the KEKB e+e" col- 
lider. Since the number of signal candidate events is compatible with that expected 
from the background, we set 90% confidence level upper limits on the branching 
fractions in the range (1.9 - 3.5) X 10-'^ for various decay modes r £ 
where I represents e oi fj,. 

Key words: TAU Lepton Flavor Violation 
PACS: 13.35.Dx, 14.60.Fg, 11.30.Hv 



1 Introduction 



In the Standard Model (SM), lepton-flavor- violating (LFV) decays of charged 
leptons arc forbidden or highly suppressed even if neutrino mixing is taken 
into account [1]. LFV is expected in many extensions of the SM such as SUSY 
models with Higgs mediation [2], right-handed neutrinos [3,4], multi-Higgs 
bosons [5], extra Z' gauge bosons [6] and i?-parity violating interactions [7]. 
Some of those models predict LFV decays of charged r leptons enhanced to 
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a level accessible at present S-factories. Observation of LFV would provide 
evidence for new physics beyond the SM. 

In this paper, we report on a search for six LFV r decay modes: t~ — > e^e^e" , 
T" —>■ e~fi'^ii'', T~ — > e'^ii~ii'', T~ — > /x~e+e~, r~ fi~^e~e~ and t~ — > 
Charge conjugate decay modes are implied throughout the paper. 
Upper limits on the branching fractions for these decays at the level (1 — 2) x 
10~^ at 90% confidence level have been set by the CLEO Collaboration using 
a data sample of 4.79 fb"^ [8]. Recently these results were improved on by the 
BaBar experiment which reported upper limits in the range (1.1 — 3.3) x 10~^ 
from a 91.5 fb^^ data sample [9]. We present here a new search based on a 
data sample of 87.1 fb^^ corresponding to (79.3 ± 1.1) million r-pairs collected 
with the Belle detector [10] at the KEKB asymmetric energy e'^e~ collider [11] 
operating at a center-of-mass energy y/s ~ 10.6 GeV. 



2 Event Selection 

The Belle detector is a general purpose detector with excellent capabilities for 
precise vertex determination and particle identification. Tracking of charged 
particles is performed by a three-layer double-sided silicon vertex detector 
(SVD) and a fifty-layer cylindrical drift chamber (CDC) located in a 1.5 T 
magnetic field. Charged hadrons are identified by means of dE/dx from the 
CDC, signal pulse-heights from aerogel Cerenkov counters (ACC), and timing 
information from time-of- flight scintillation counters (TOP). Energies of pho- 
tons are measured using a Csl (Tl) electromagnetic calorimeter (ECL). Muons 
are detected by fourteen layers of resistive plate counters interleaved with iron 
plates (KLM). 

We search for t~^t~ events in which one r decays into three charged leptons 
(3-prong) and the other r decays into one charged and any number of neutrals 
(1-prong), which has a branching fraction of (85.35±0.07)% [12]. We start by 
requiring signal candidate events to have four charged tracks with zero net 
charge and any number of photons. Each charged track must have transverse 
momentum pt > 0.1 GeV/c and be within the polar angle range 25° < 9 < 
140°. For each charged track, the distance of closest approach to the interaction 
point (IP) is required to be within ±1 cm transversely and ±3 cm along 
the beam. Photon candidates are selected from neutral ECL clusters with an 
energy -Eduster > 0.1 GeV, where neutral ECL clusters must be separated by 
at least 30 cm from a projection point of any charged track in the ECL. The 
tracks and photons in an event are divided into two hemispheres in the e+e~ 
center-of-mass system (CMS) by the plane perpendicular to the thrust axis 
calculated from the momenta of all charged tracks and photons in the event. 
Signal candidates have a 1-prong vs 3-prong topology, i.e. three charged tracks 



4 



are required in one hemisphere, one charged track in the other. We define the 
former hemisphere as the signal side and the latter as the tag side. The number 
of photon candidates on the signal side, n^, should be less than or equal to 
two, to allow for photons from initial state radiation or photon radiation in 
the detector by electron tracks. 

Electrons are identified by means of an electron likelihood function (Ce), that 
includes the information on the dE/dx measurement in the CDC and the ra- 
tio of the cluster energy in the ECL to the track momentum measured in the 
CDC [13]. For electrons, we require >Ce > 0.1 and the laboratory momentum 
to be greater than 0.3 GcV/c. In order to correct for the energy loss from 
bremsstrahlung in the detector material, the momentum of an electron can- 
didate is recalculated by adding the momentum of radiated photon clusters 
when an ECL cluster with energy less than 1.0 GeV is detected within a cone 
angle of 10° around the flight direction of the electron candidate track. The 
muon likelihood function is evaluated from two variables - the difference 
between the range calculated from the momentum of a particle and the range 
measured in the KLM as well as the value of the KLM hits with respect 
to the extrapolated track [14]. For muons, is required to be larger than 
0.1 and its momentum should be larger than 0.6 GeV/c. The requirement of 
lepton identification for all three charged tracks on the signal side leads to 
significant background reduction by a factor in the range from 10^ to 10^ for 
each decay mode with a 40% loss of signal. 

For the Monte Carlo (MC) simulation of the signal, one r lepton is assumed to 
decay into the signal LFV modes with a uniform phase space distribution in the 
r rest frame and the other r decays generically using the KORALB/TAUOLA 
program [15]. 

One important source of background remaining after requiring the event topol- 
ogy and lepton identification on the signal side is due to radiative Bhabha 
events with a converted photon. This background is efficiently reduced by re- 
quiring that the invariant mass of any two oppositely charged particles be 
greater than 0.2 GeV/c^, assuming the electron mass for each particle. The 
remaining background after removing photon conversions comes from two- 
photon processes, r-pair events with generic decays into three charged hadrons, 
e+e" qq continuum and BB events, where hadrons are misidentified as lep- 
tons on the signal side. 

For signal r-pair events there is a missing momentum due to neutrino emission 
from the r on the tag side. Figure 1 (a) shows the total transverse momentum 

p^, the transverse component of the sum of momentum vectors for the four 
charged tracks in the CMS. To suppress events from two-photon processes, 
is required to be larger than 2.0 GeV/c. We calculate the missing momentum 
by subtracting the momentum of all charged tracks and photons from the beam 
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Exp.data I ! signal MC two-photon ^ I Bhabha and ^-pair others 




Fig. 1. Kinematical distributions used in the event selection after the event topology 
and lepton identification requirements: (a) the total transverse momentum p^, (b) 
the polar angle of the missing momentum ^miss, (c) the opening angle between the 
momentum of a 1-prong track and the missing momentum ^fp-migs, (d) the momen- 
tum of a 1-prong track pjp. The distributions for experimental data, signal MC and 
background MC are indicated by closed circles, the dashed and solid histograms, 
respectively. The areas of the signal MC are normalized assuming a branching frac- 
tion of 1.3 X 10^^ while the background MC is normalized to the data luminosity. 
Background MC from two-photon, Bhabha and /i-pair, and other processes (generic 
rr, e^e~ — > qq continuum and BB) are shown by the open, shaded and dark his- 
tograms, respectively. The distributions are shown for the r~ e~e~^e~ mode. 
Although the size of each background source depends on the decay mode, the shape 
of each background is common to all. The shape of signal distributions is also similar 
for the six decay modes. 

momentum. In order to suppress radiative Bhabha and two-photon events, the 
polar angle of the missing momentum in the laboratory frame ^miss, niust be 
between 25° and 140°, as shown in Fig. 1 (b). The missing momentum due to 
neutrinos from the r on the tag side tends to lie in the same hemisphere as the 
1-prong track for signal events. The opening angle between the 1-prong track 
and missing momentum in the CMS, ^ip_miss; is required to be less than 90°, 
as shown in Fig. 1 (c). Finally, the momentum of the 1-prong tag side track, 
pip, must be less than 3 GeV/c. Since the tag side r decays with neutrino(s) 
or 7r° emission, the momentum of the 1-prong track is much smaller than 
the T momentum. This requirement suppresses most of the Bhabha and /x- 
pair backgrounds, as shown in Fig. 1 (d). After all these requirements, the 
background is reduced by a factor of order 10^ with a 10% efficiency for the 
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signal. 



From the total CMS energy, E^^, and invariant mass, M^g, of the three signal 
side leptons, we compute the quantities: AE* = E^^^ — -Ebeam ^'^^ = 
M-^i — Mr, where -E'^eam is the CMS beam energy and M,- is the r mass. In the 
AE*-AM plane, the neutrinolcss r decay events are expected to be distributed 
close to the origin. The AE* and AM distributions for each decay mode are 
shown in Figs. 2 and 3, respectively together with expectations based on signal 
MC. It is seen that each peak has a tail on the lower side that is due to initial 
state radiation or bremsstrahlung of a charged particle interacting with the 
detector material. Since electrons radiate more than muons, the shape of the 
peak depends on decay mode. 

The signal region for each decay mode is given in Table 1, and is illustrated 
as the region between two dashed lines in Figs. 2 and 3. Each signal region is 
defined to contain 90% of the signal MC events plotted in the figures. Figure 4 
shows the AE* vs AM plot for the data. In the fourth and fifth columns 
of Table 1, we compare the number of events found for the data and the 
normalized MC background in the plotted region of —0.68 GeV < AE* < 
0.32 GeV and -0.12 GeV/c^ < AM < 0.12 GeV/c^ shown in Fig. 4. 



3 Results 

Efficiencies for r — > 3£ decays with a uniform phase space distribution vary 
from 9.2% to 9.5% and are listed in the second column of Table 2. The actual 
decay angle distribution, however, will depend on the LEV interaction and 
include spin correlations between the tag side and signal side r [16]. In order 
to evaluate the maximum possible effect of such correlations, we examine 
V — A and V + A interactions using the formulae given in [16]. The relative 
differences in the efficiencies (Ae/e) from a uniform distribution are found to 
be -3.8%, -8.7%, -1.1%, +0.8%, -12.6% and -5.6% for r" ^ e-e+e", r" ^ 
e~ li~ , e~^fi~fi~, T~ — >• /x~e~e"'", r~ — >• ii^e~e~ and r~ — >• ii~ ii^ ii~ 

decay mode, respectively. Requirements on the number of CDC tracks and the 
energy of ECL clusters that are used to detect r-pair events also constitute part 
of the trigger logic. The impact of the trigger on the efficiency is investigated 
by applying a trigger simulation to the signal MC. The changes in the detection 
efficiencies are found to be only a few percent, because the selection criteria 
applied in this analysis are much more restrictive than the trigger conditions. 

As shown in Fig. 4, the background level in and around the signal region is 
very low. From the background MC, we find that the remaining events are due 
to the low-multiplicity e+e~ — > qq continuum or BB events where final-state 
hadrons are misidentified as leptons as well as a few events from generic r-pair 
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Fig. 2. AE* distributions of the events in the AM signal region. Open histograms 
are the signal MC sample assuming the branching fraction of 1.3 x 10~^, while 
experimental data is plotted by closed circles. The dashed lines indicate the AE* 
signal region. 
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Fig. 3. AM distributions of the events in the AE* signal region. Open histograms 
are the signal MC sample assuming the branching fraction of 1.3 x 10~^, while 
experimental data is plotted by closed circles. The dashed lines indicate the AM 
signal region. 

decay. As seen from Table 1, the numbers of events found in the Fig. 4 area 
are consistent with the numbers expected from the normalized background 
MC. 

To evaluate the background b in the signal region, we assume a uniform back- 
ground distribution along the AM axis in Fig. 4. With looser selection cri- 
teria, we find the AM distribution is uniform. We estimate the number of 
background events in the signal region from the number of events observed in 
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Fig. 4. Ai?* vs AM plots for the experimental data. The charge conjugate decay 
mode is also included. The dashed and dotted boxes indicate the signal and AM 
side-band regions, respectively. 

Table 1 

Definition of the signal regions for each decay mode and the number of events in 
Fig. 4 for the data and the normalized background MC. 

Mode Signal region Number of events 

in the Fig. 4 area 











A^;* (GeV) 


AM (GeV/c2) 


Data 


MC 




e' 




e 


-0.36 < A^;* < 0.04 


-0.032 < AM < 0.010 


1 


u -0 


- 


-> e~ 






-0.32 < A^;* < 0.03 


-0.017 < AM < 0.010 


18 


5 +17 


- 








-0.32 < AE* < 0.03 


-0.017 < AM < 0.010 





9 +3 
^ -2 




fJ-' 




-e+ 


-0.33 < A^;* < 0.04 


-0.025 < AM < 0.010 


2 


U -0 


r~ 


fJ-' 






-0.33 < A^;* < 0.04 


-0.025 < AM < 0.010 





+3 
u -0 


- 


f^" 






-0.28 < AE* < 0.03 


-0.010 < AM < 0.010 


5 


8±4 



the dotted box, the AM side-band regions in Fig. 4. The systematic uncer- 
tainty for this method is estimated by comparing the observed and estimated 
number of events in the region outside of the dotted box in Fig. 4. The fourth 
column of Table 2 shows the number of estimated background events and its 
error, including both systematic and statistical uncertainties. 
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Table 2 

Summary of detection efficiency, number of observed events, background expecta- 
tion, So and 90% C.L. upper limits on the branching fractions. 



Decay 


Detection 


Number of 


Expected 




Upper limit 


mode 


efficiency 


events observed 


background 


•So 


of ^ 




(e), % 








xlO"^ 


T~ — ^ e~e'^e'' 


9.2 ±0.6 


1 


< 0.2 


4.36 


3.5 


T" — e~ fjL" 


9.2 ± 1.4 





0.1±0.1 


2.54 


2.0 


T~ — > ii~ jjL" 


9.2 ± 1.1 





< 0.3 


2.55 


2.0 


T~ — >■ //~e~e+ 


9.4 ±0.8 





0.2±0.2 


2.49 


1.9 




9.5 ± 1.4 





< 0.2 


2.55 


2.0 


T~ — > pL~ pL^ ll~ 


9.0 ± 1.6 





0.1±0.1 


2.51 


2.0 



In the signal regions for the six decay modes considered, one candidate is 
observed for the t~ e~e^e~ mode while no candidates are found for the 
other modes. The numbers of events in the signal regions are consistent with 
the background expectations. 



We determine the upper limit Sq on the number of signal events at 90% CL 
using the prescription of Feldman and Cousins [17]. To include in this limit 

the uncertainty in the detection efficiency e, we increase sq according to the 
prescription of Cousins and Highland [18]. The main systematic uncertain- 
ties in the detection efficiencies come from tracking (1.0% per track), electron 
identification (1.1% per electron), muon identification (5.4% per muon), trig- 
ger efficiency (1.4%), statistics of signal MC (1.0%) and uncertainty of the 
decay angular distribution (0.8 — 12.6%). The total uncertainty of the detec- 
tion efficiency is 6.1% for r~ e~e+e^, 15.1% for e^fi^fj,^, 12.4% for 
T~ —>■ e~^fi~fi~, 8.4% for r~ — >• /i~e~e+, 15.1% for r~ — >• ii'^e~e~ and 17.5% for 
T~ —>■ ii~ mode. The uncertainty in the number of r-pair events comes 
from the luminosity measurement (1.4%). For our calculation of sq, we take 
the background (and its uncertainty) to be zero. This results in conservative 
upper limits. Upper limits for branching fractions B are calculated for each 
decay mode as follows: B{t~ — > i~£^i~) < — , where N^-^- is the 

total number of the r-pairs produced, and Bi is the inclusive 1-prong branch- 
ing fraction of the r. The values of Sq used and the resulting upper limits for 
the branching fractions are summarized in Table 2. 
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4 Summciry 



We have searched for lepton-flavor-violating decays t~ — > t't^t' using an 
87.1 fb^^ data sample. No evidence for any of these decay modes is observed 
and upper Umits for the branching fractions are obtained in the range (1.9 — 
3.5) X 10^^ for I'l'^t' modes; these are approximately one order of 

magnitude more restrictive than the hmits previously obtained by CLEO [8] 
and comparable to the recent results from BaBar [9]. 
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